Strong-coupling branching of FQHL edges 



^1" 

O 

o 

(N 

C 

in 






c3 



13 

o 
o 



> 

oo 

O 

O 

-i— > 

c3 



i 

C 

o 
o 



X 



V.V. Ponomarenko and D.V. Averin 
Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, NY 11794, USA 

(February 2, 2008) 

We have developed a theory of quasiparticle backscattering in a system of point contacts formed 
between single-mode edges of several Fractional Quantum Hall Liquids (FQHLs) with in general 
different filling factors i/j and one common single-mode edge j/q of another FQHL. In the strong- 
tunneling limit, the model of quasiparticle backscattering is obtained by the duality transformation 
of the electron tunneling model. The new physics introduced by the multi-point-contact geometry 
of the system is coherent splitting of backscattered quasiparticles at the point contacts in the course 
of propagation along the common edge v®. The "branching ratios" characterizing the splitting 
determine the charge and exchange statistics of the edge quasiparticles that can be different from 
those of Laughlin's quasiparticles in the bulk of FQHLs. Accounting for the edge statistics is 
essential for the system of more than one point contact and requires the proper description of the 
flux attachement to tunneling electrons. 



I. INTRODUCTION 

Electron transport properties of the Fractional Quan- 
tum Hall Liquids (FQHLs) continue to attract consid- 
erable interest motivated by the unusual properties of 
FQHL excitations which are characterized by fractional 
charge and non-trivial exchange statistics. Since exci- 
tations in the bulk of a typical FQHL are suppressed 
by the energy gap, the low-energy transport properties 
are determined by the gapless modes at the edges of the 
liquid 1 . In the simplest case of a FQHL with fillings fac- 
tor v = 1/odd, its sharp edge should support one bosonic 
mode that is described as a single branch of chiral Lut- 
tinger liquid 2 . Although the edge states are gapless and 
in principle can carry excitations of arbitrary charge, in 
the case of tunneling between the edges of FQHLs with 
equal filling factors, strong tunneling produces excita- 
tions which coincide 3 with Laughlin's quasiparticles in 
the bulk of the liquid. Experiments on resonant tunnel- 
ing through a quantum antidot and on the shot noise 
in quasiparticle backscattering in a single point contact 5 
measure directly the fractional charge of these excita- 
tions. Edge-state tunneling should also make possible 
the measurements of quasiparticle exchange statistics 6 
which can be used in the development of FQHE qubits 7 
for solid-state quantum computation. 

In the situation of tunneling between the edges of 
FQHLs with different filling factors, strong tunneling 
should produce quasiparticles that are different from bulk 
quasiparticles 8 . The charge of such "contact" quasipar- 
ticles coincides with the dc conductance of the point 
contact (if the two are measured in units of electron 
charge e and the free-electron conductance e 2 /h, respec- 
tively.) Theory up to now could predict only the prop- 
erties of one point contact between different edges. The 
aim of this work is to develop a theory of strong tun- 
neling between FQHL edges with different filling factors 
in a junction with more than one point contact. Such a 
model of multi-point-contact junction was introduced 9 ' 10 



to describe experiments 11 on tunneling between FQHL 
edge and external Fermi-liquid reservoir. Understand- 
ing of the strong-tunneling limit of this model is impor- 
tant for the description of the reservoir-edge equilibration 
that is responsible for correct quantization of the two- 
terminal FQHL conductance. 12 Multi-point-contact junc- 
tions with a well-controlled geometry similar to the one 
considered in this work are also studied experimentally. 13 

The main technical obstacle to the description of 
strong multi-point-contact tunneling between different 
edges is the fact that the bosonic fields in the tunnel- 
ing operators of different contacts do not commute with 
each other and, at first sight, can not be localized simul- 
taneously at the minima of large tunnel potential, as can 
be done for one point contact. This problem is resolved, 
however, if the bosonic fields are modified 12 by the proper 
choice of the statistical phase of the tunneling electrons 
[see Eqs. (11) and (12) below]. The statistical phase pre- 
serves the Fermi statistics of electrons but accounts for 
the change of the number of flux quanta attached to them 
in the FQHLs 14 in the process of tunneling between the 
liquids with different filling factors. With the flux at- 
tachement taken into account, the total bosonic tunnel- 
ing fields commute, and one can build the strong-coupling 
description of the multi-point-contact junction in close 
analogy to the case of one point contact. Here, we use 
this approach to develop complete description of strong 
tunneling in the multi-point-contact junction extending 
previous results 12 to quasiparticle backscattering. 

Important elements of our approach can be summa- 
rized as follows. Specific junction model considered in 
this work is characterized by the existence of one edge i>q 
common to all point contacts (Fig. 1). The quasiparticle 
backscattering is produced by finite reflection coefficients 
of the contacts, and is described as instanton tunneling 
between the infinite set of the ground states of the orig- 
inal electron tunneling model that arc degenerate in the 
absence of backscattering. Expansion of the junction dy- 
namics in terms of instantons gives the model of quasi- 



particle tunneling which is dual to the electron tunneling 
model. Duality transformation relating the two models 
produces the quasiparticle exchange statistics from the 
electron statistical phases ascribed through the flux at- 
tachemcnt. Correct description of quasiparticle statistics 
in our junction geometry enables one to combine sequen- 
tially the edge-state transformations at different point 
contacts in such a way that the incoming edges are split 
coherently at each point contact in the course of propa- 
gation along the common edge vq. This process can be 
described with the "branching ratios" which define how 
the incoming currents (or charges) are split between the 
outgoing edges. Quasiparticles produced by backscatter- 
ing at one point contact are split similarly at all point 
contacts downstream the edge vq from the "contact of 
origin" , and the edge-state branching ratios determine 
the charges generated in each edge by backscattering. As 
usual, these charges should manifest themselves through 
the intensity of the shot noise of outgoing currents. The 
charge splitting also provides natural interpretation of 
the quasiparticle exchange statistics obtained by the du- 
ality transformation, so that both the charge and statis- 
tics of backscattered quasiparticles can be viewed as be- 
ing defined by the strong-coupling edge-state branching. 
The paper is organized as follows. Section II defines 
the electron tunneling model considered in this work. 
Section III describes the bosonization procedure for the 
Klein factors of electron tunneling operators that im- 
plements the flux attachment. Section IV explains the 
instanton transformation, formulates the dual model of 
quasiparticle tunneling, and calculates the dc current and 
shot noise in outgoing edges of the junction. In Section 
V we apply the results of Sec. IV to the situation that 
corresponds to the junction between a Fermi- liquid reser- 
voir and FQHL edge. Quasiparticle backscattering in 
this case gives rise to corrections to the Ohmic behavior 
of the junction: non-linear current-voltage characteristic 
and shot noise, which limit the reservoir-edge equilibra- 
tion obtained in the absence of backscattering. In the 
Appendix, we summarize for convenience the known re- 
sults for strong tunneling in one point contact that are 
used in the main text. 



II. MULTI-EDGE MODEL OF WEAK ELECTRON 
TUNNELING 

The model we consider consists of n edges of different 
FQHLs of the filling factors i/j, j = 1 ...n, which form 
n tunnel point contacts with an edge of another FQHL 
of the filling factor vq (see Fig. 1). All liquids are as- 
sumed to correspond to simple Laughlin states with the 
filling factors v\ = 1/odd, 1 = 0... n, which in general can 
all be different. Since bulk excitations inside the liquids 
have energy gap, at energies below some common cut-off 
energy D, only the edges support excitations that prop- 
agate along them. In the weak tunneling regime, when 
each of the jth edges is well separated from the edge, 



the tunnel currents between them should be carried by 
individual electrons. We assume that the tunnel contact 
formed by the jth edge is located at the point Xj along 
the edge and its width is much smaller than the mag- 
netic length. Tunneling in such a system of contacts is 
described by the Lagrangian: 
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j'=i 



[UjtpQ (xj , t)ipj (xj , t) + h.c. 



(1) 



where the tunnel amplitudes Uj are taken to be real and 
positive, and ipi is the electron operator of the Ith edge. 
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FIG. 1. Diagram of the model of a tunnel junction consid- 
ered in this work: n edges of FQHLs with in general different 
filing factors Vj tunnel into one FQHL edge with the filing fac- 
tor VQ at points Xj along this edge. Uj are electron tunneling 
amplitudes. Each edge is assumed to support one branch of 
bosonic excitations with arrows indicating direction of prop- 
agation of these excitations. 

If we assume that all edges are sharp enough to enable 
us to use simple bosonic single-mode description of the 
edge-state dynamics introduced by Wen 2 , the electron 
operator tpi can be expressed as 



ipl = {2ira)- 1/2 ^e iMx ' t)/ ^ U ' 



(2) 



Here <j>i is the bosonic field of the edge I; Majorana 
fermions £/ account for mutual statistics of electrons in 
different edges, and a common factor 1/a = D/v denotes 
momentum cut-off of the edge excitations. Since the spa- 
tial dynamics of the edges j = 1 ... n does not affect the 
tunnel currents, propagation velocities of these edges are 
irrelevant, and we take them to be equal to the velocity 
v of the excitations of the edge v$. The bosonic fields 4n 
are normalized in such a way that their free dynamics is 
governed by the standard Lagrangian density 2 
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Y-d x (d t +vd x ), 



which corresponds to the Hamiltonian 

Ho = ^ J dx(d x Mx)) 2 
and the commutation relations 



(4) 



[4>i (x),</> p (0)]=iw 5i p sgn(x) . 



(5) 



In the following, we will need free correlator of the 
fields 4>i ordered in imaginary time r: 

Gi p (x,t) =6i p g(x,T), g(x,r) = (T t {4>i(x,t)<Pi (0,0)}) . 

(6) 

To find this correlator, we start with the retarded Green's 
function in real time t: 



g R (x,t) = -*e(t)<[Ma,i),0,(O,O)]> , 



(7) 



where 0(t) is the step function. Lagrangian (3) implies 
that dynamics of g R (x, t) is governed by the equation: 

d x (d t + vd x )g R {x,t) = 27rS(x)5(t) . 

Solution of this equation that satisfies the symmetry of 
the Eq. (7) (the commutator part of Eq. (7) should be an- 
tisymmetric with respect to simultaneous change of sign 
of x,t) is: 

g R (x, 1) = 7r6(t)sgn(x - vt) . 

The standard steps: Fourier transformation of g (x,t), 
g R (x,oj) = J dte luJt g R (x,t), and analytical continuation 
to the upper half-plane of frequency, lo — > ilo, give the 
Fourier transform of the imaginary-time correlator for 
positive frequencies lo: g(x,Lo) = — g R (x, lo)]^^^. 

Similar calculations for the advanced Green's function 

g A (x, t) = i9(-t)([4>i(x, f ). M) = -7re(-t)sgn(a; - vt) , 

give g(x, lo) for negative frequencies, and it can finally be 
written for all Matsubara frequencies as 



l(x,u) = ^sgn(x) (- 1 - + Q(wx)e- ux ' v \ 



(8) 



The first part of this correlator does not depend on |x| 
and characterizes singular behavior of the correlator at 
short times, g(x,r) <~ i7rsgn(xr)/2, which is responsible 
for the intra-edge statistics of the electronic operators 
ipl. With the adopted normalization of <j>i, the operator 
of charge density of the Ith edge at point x is 

pl{x,r) = ( v fv~i/2ir)d x ({>i(x,T) , 

and the current in the edge is related to the charge den- 
sity as ji = vpi. 

As in the case of a single point-contact, the effec- 
tive electron tunneling amplitudes are renormalizcd from 
their original values Uj in Eq. (1) by fluctuations of the 
fields <f>i. The effective tunneling amplitudes scale with 
energy in such a way that the weak-tunneling limit is al- 
ways stable at sufficiently low energy. The main focus 
of this work is on the strong-tunneling regime realized at 
large temperatures or bias voltages Vj across the point 
contacts which make the effective tunneling amplitudes 



large. The description of the strong-tunneling limit pre- 
sented below is based on the instanton expansion around 
the saddle-point solution 12 . As discussed in the Intro- 
duction, this solution requires appropriate bosonization 
of the Klein factors of the electron tunneling operators 
that takes into account the flux attachcmcnt. 



III. BOSONIZATION OF THE KLEIN FACTORS 

Substitution of the bosonic representation (2) of the 
electron operators tpi into Eq. (1) transforms the tunnel- 
ing Lagrangian into 



t^tunn — / , I -. 



Fj exp {iXj(fj{t)} + h.c. 



(9) 



Here 



\jipj(t) 



and the factors 



<f>o(xj,i) <j>j(xj,t) 



"I 1/2 



A; = 



V Q + Vj 



V Vj 



(10) 



are chosen in such a way that the normalization of the 
bosonic operator (fj coincides with the normalization of 
the bosonic fields 4>i used before, so that the imaginary- 
time correlator of tpj is given by the same Eq. (8) with 
x = 0: g(0,uj) =n/\u>\. 

In the Lagrangian (9), Fj represents the anticommut- 
ing statistical Klein factors, Fj — £,o£,j, which account for 
mutual statistics of electrons in different edges. These 
factors play an important role in the strong multi-edge 
tunneling considered in this work. A convenient way of 
taking them into account is to express them through the 
zero-energy bosonic fields r\j : 12 ' 15 



Fj = e^ , 



[Vi,Vj] =mi%3 



(11) 



where 7^ are odd integers. Different choices of 7^ corre- 
spond to different branches of the phase of the Fermionic 
statistical factor —1 arising from interchange of Fi and 
Fj . Under the condition that the point contacts are well- 
separated, \xi — Xj\ 3> a, minimization of energy of the 
strong-coupling ground state requires that the statistical 
phases are taken as: 12 



lij =sgn(i -j){l-Sij) 



1 



(12) 



This equation assumes that the point contacts are num- 



bered as in Fig. 1: 



< Xj for 1 < % < j < n, and 



the coordinate x increases in the direction of the edge 
propagation. 

The steps similar to those that lead to Eq. (8) give 
then the Fourier transform of the r-ordered correlator of 



the fields 



Vj- 



(ViVj) =sgn(t-j)(l -6 i:j ) 



uv n 



(13) 



Substitution of the bosonizcd form (11) of the Klein fac- 
tors into Eq. (9) turns the tunneling Lagrangian into a 
function of n bosonic variables 



<bj = Xjipj 



Vj 



(14) 



and one can reduce the kinetic Lagrangian (3) to a sim- 
plified form by integrating out all the fields except <F,. 
The resulting Gaussian action for the fields $j can be ob- 
tained directly by noticing that it should be determined 
by the matrix of correlators -Ky(u>) = ($»(— u))$j(uj)) as 

So = U&- 1 $= ^^(-c^V^-M, (15) 



where lo = 2i:mT denotes the Matsubara frequencies 
with m — ± 1,±2, ... and temperature T. The ma- 
trix of correlators K^ can be found from the correlator 
(8) of the tunneling fields ipj (determined by the original 
kinetic Lagrangian (3)) and the correlators (13) of the 
statistical fields. Combining the two we get 






2 



n 



9(sgn(i-j)w) _ 

e 



Wti. 



'(!■ 



Sij) 



(16) 



where tij — U — tj and tj = Xj/v. Important feature of 
the terms in Eq. (16) that are non-diagonal in i, j is that 
they do not contain the singular statistical parts. The 
singular part of the correlator (8) of the tunneling fields 
ifj is cancelled by the statistical fields r]j, if the phases 
jij (11) are chosen appropriately as in Eq. (12). 



IV. STRONG-TUNNELING LIMIT 

A. Instanton expansion and duality transformation 

The tunneling action (9) expressed through the bosonic 
fields $j (14) has an infinite set of minima at $j = 
27r x integer, j = 1 -f- n, that are degenerate in en- 
ergy In the strong-tunneling limit, when the electron 
tunneling amplitudes Uj are large, these minima are well- 
separated in a sense that the amplitude of 3>j -tunneling 
between them is small. Such rare tunneling processes 
correspond physically to finite backscattering at each of 
the point contacts j and are described by the instanton 
tunneling solutions which on the long-time scale behave 
as 2ireijO(T — ry), where index I counts different instan- 
ton tunneling events, eij = ±1, and tij are the times of 
tunneling of the <&j component. On the short time scales 
on the order of 1/D, the variation of <J?j is smoothed 
around ry, with the exact behavior dependent on the 
form of the energy cut-off. For some special form of the 
cut-off, 17 the shape of the instantons for short times can 



be found from the equations of motion. It minimizes the 
instanton contribution to the action in the absence of 
the long-time interactions and determines the instanton 
tunneling amplitudes. The amplitudes are, in addition, 
modified by quantum corrections and we take them as 
some unspecified parameters Wj/(2ira). 

The low-energy part of the action determines the 
instanton-instanton interaction and can be found by sub- 
stitution of the long-time asymptotic form of $j(r), i.e., 
$j(t) = 2tt x integer + ^ ; ey27r6(r — Ty) into the (In- 
action (15). Indeed, the constant part of $j(t) is irrele- 
vant, and as we will see later, only the trajectories with 
the same initial and final values of <&j(t) for each con- 
tact j, ^ ; eij — 0, arc important. The Fourier transform 

$j(w) = VT jl /T dTe lu)T ^ 1 {T) of such "neutral" trajec- 



tories is: 



*» 






eij- 



2ms/T 



(17) 



and combining this equation with Eq. (15) we see that the 
low-energy part of the action has the form of the pair- 
wise instanton interaction with each pair of instantons 
with indices Ij and ki contributing the term 



S int (Ti,T k ) = e l3 e kl TY, (^) K^{uj)e^ Tl - Tk \ 



(18) 



It is convenient to formulate the instanton dynamics 
in terms of the fields Qj, j — 1 ... n, that are defined as 
dual to <I>j, i.e., satisfy the commutation relations: 



[%^ 



2Tri5j 



(19) 



In terms of these fields, each instanton tunneling is gen- 
erated by the operator 



2^exp{ie y e j -('n)} 



(20) 



which shifts &j by 2ireij . The low-energy part of the ac- 
tion containing instanton-instanton interaction (18) can 
be understood then as arising from Gaussian fluctuations 
of the dual fields: 

e -S int _ e -e, J e ki (e l (r k )e j (n))/2 
and Eq. (18) gives the correlators of these fields: 
(Q i (-u)& j (u;)) = (2Tr/u,fK- j 1 (u;). 

To find these correlators explicitly, we need to invert 
the matrix K defined in Eq. (16). This task is not diffi- 
cult, since the matrix has a triangular form: for negative 
frequencies, K and K~ x are upper-triangular, while for 
positive frequencies they are lower-triangular. In both 
situations, one can find K^ 1 by writing equations for the 
matrix elements Kj~- as recurrence relations in terms of 



the "distance" from the diagonal and solving them start- 
ing with the diagonal elements. For negative frequencies, 
the result is: 



2tt 



^Sij - e(sgn(i-j) W )B«e-l wt «l(l - %) 



(21) 



The representation (24) of Qj allows us to write the in- 
stanton tunneling Lagrangian in the form analogous to 
Eq. (9): 

Ctunn = Y^-Fj exp{*-?:^(*)} + h.c] . (26) 



i=i 



2ira 



A, 



where the Klein factors of the dual model are: 



Bi 



n 



^oAf A2 . , 

4 J i<k<j 



VoH 



(22) 



F,- 



=*% 



[VuVj] 



-i7rsgn(i — j)(l 



S ij) B i. 



(27) 



For positive frequencies, Eq. (21) remains valid if we de- 
fine the matrix Bij for i > j (below the diagonal) by 
the condition B^j = Bji. Qualitatively, as will be dis- 
cussed in more details below, the dual fields Qj describe 
quasiparticles backscattered at the point contacts, and 
the matrix B^ (22) determines the exchange statistics of 
these quasiparticles. 

Similarly to the case of electron tunneling model where 
the correlators (16) give the kinetic part (15) of the ac- 
tion, the correlators (21) determine the kinetic part of 
the action of the dual model of quasiparticlc tunneling: 



S = ^Q{w/2ir) 2 KQ . 



(23) 



One more remark is that the 9 j -correlators (21) show 
that the condition of neutrality of instanton tunneling 
trajectories, ^ ; eij = 0, that was used above, is satisfied 
at each point contact j. Indeed, the form (21) of the 
correlators implies that interaction between instantons 
diverges at low energies and only the trajectories that 
satisfy the neutrality condition for each j have finite ac- 
tion and contribute to the evolution of the system. 



B. Dual chiral fields and Klein factors 

One can see directly that the expansion of the system 
propagator in the number of instanton tunneling events 
can be generated as the expansion in the instanton tun- 
neling Lagrangian composed of the tunneling operators 
(20). This means that the dual model of quasiparticle 
tunneling can be formulated in complete analogy to the 
direct model of electron tunneling. Comparing the corre- 
lators (21) with those from Eq. (16), one can notice that 
similarly to the $ fields (14), Qj can be represented as 
the sum 






m 



(24) 



of the fields with, respectively, chiral and pure statistical 
correlators. Equation (21) implies then that the chiral 
correlator of the fields Qj for i ^ j is: 



(6i6j) = -{\i\j/A)Bijg{xi - Xj,u>) . 



(25) 



Before discussing the physical consequences of the dual 
model of quasiparticle backscattering, we present the 
derivation of this model that is less rigorous than the one 
that uses instanton expansion, but more direct. The ap- 
proach is based on independent application of the known 
strong-tunneling solution for one point contact 8 to indi- 
vidual contacts of our multi-point-contact junction and 
matching the obtained solutions at successive contacts 
along the edge vq common to all of them (see Fig. 1). 
Since the bosonic fields at different point contacts inter- 
act strongly at low frequencies, the possibility to simply 
match the solutions at successive contacts is by no means 
trivial and represents an important assumption. The in- 
stanton calculation described in the previous subsection 
can be viewed as the proof that this assumption is in- 
deed valid provided one chooses correctly the statistical 
phases of electrons in different contacts, which after du- 
ality transformation determine the statistical phases of 
backscattered quasiparticles. 

In more details, consider the jth point contact of 
our multi-contact model. Strong-tunneling solution (re- 
viewed in the Appendix) can be described as application 
of the Dirichlct boundary condition to the tunneling field 
<fij{x) which localizes this field at the minimum of the 
tunnel potential at the point of contact j. "Unfolded" 
form of this condition 18 implies free propagation of the 
two fields constructed from the original bosonic modes 
4>j(x) and (j>o(x) of the edges forming this contact. One 
is the field 



9j(x) = — sgn(x — Xj)(pj(x) 



(28) 



dual to the tunneling field (pj(x), while the other, <fj{x) 
is the combination of <f>j (x) and (f>o (x) orthogonal to the 
tunneling field, i.e., 



1/00 <t>3 \ - 1/00 



(29) 



(We use the same notation for the dual field Qj defined 
by Eq. (28) as in Eq. (24) in anticipation of the fact, 
proven later, that these two fields indeed coincide.) Free 
evolution of the field Qj(x) (28) is equivalent to the state- 
ment that the incoming tunneling field changes sign at 
the contact x = Xj . Combining this change of sign with 
Eq. (29), we see that the transformation of the incom- 
ing fields (f>o (x) and 4>j{x) at x < Xj into outgoing fields 



at x > Xj (Fig. 2a) can be described 8 with the 2-by-2 
matrix P«': 



H^M^pW^H^W)] 



with the matrix elements 



50) 
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0) _ ^o - ^ 60) _ 60') _ 2 V 77 " 



i/ + I/j 



P, 
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p 



J'O 



(out) 



I/Q + Z/j 



(30) 
(31) 




FIG. 2. Diagram of the strong-coupling scattering of the 
edge states in the multi-point-contact junction considered in 
this work: (a) transformation of the bosonic edge-state fields 
(30) in one individual point contact j; and (b) transformation 
of the fields in the n-point junction as a whole: the field <f>o 
undergoes successive transformations in the course of propa- 
gation along the edge Vo. 

A natural way to combine transformations (30) at the 

neighboring contacts is to assume that the field <fr OM 
coming out of each contact along the edge i/ serves as 
the incoming field </> for the next point contact down- 
stream this edge (Fig. 2b). Under this assumption (va- 
lidity of which implies the non-trivial mutual statistics 
of backscattered quasiparticles), description of the over- 
all transformation of the incoming bosonic modes in our 
multi-point-contact junction can be obtained by extend- 
ing each of the matrices PO') trivially to the n+1 by n+ 1 
matrix: the fields (pi with i ^ j arc constant in pv>, and 
multiplying this matrices. The full transformation of the 
incoming fields cjy- 17 ^ — (</>o, (pi, ... , 4> n ) T into the outgoing 
fields (see Fig. 2b) is described then as: 

^(out) = §<p( in \ 5 = P (n) P (n_1) ■ ... -P (1) . (32) 

Such a combination of successive transformations of 
the fields at the point contacts together with the free 
propagation of the field </>o between the contacts implies 
that the correlator of the fields Qj(x) defined by Eqs. (28) 
and (29) at different contacts is: 



i(xi)9j(xj)) = - 



1 



voXiXj 



){Xi 



"J! 1 



' n 

i<k<j 



P 



(k) 



00 



Equations (10), (22), and (31) show that this correla- 
tor coincides with the correlator (25). According to 



Eq. (28), the correlator (9j6j) at the same contact is 
(OjOj) = {tpjifij) — t/|u>| and also agrees with Eqs. (24) 
and (21). This means that the chiral fields 9j intro- 
duced in Eqs. (24) and (26) as the fields describing the 
quasiparticles backscattered at the point contacts coin- 
cide with the dual fields 9j(xj) from Eq. (28) matched 
directly between the contacts. This coincidence is non- 
trivial and is only correct if the mutual statistics of quasi- 
particles backscattered at different point contacts is ac- 
counted for by inclusion of the Klein factors Fj (27) . The 
statistics, however, does not affect the scaling dimensions 
of the quasiparticlc tunneling operators, which have the 
same values 2/X 2 as for one independent point contact 19 . 
This means that the dual model of strong tunneling in 
the multi-point-contact junction constructed in this work 
is stable for large temperatures and/or voltages applied 
to all contacts, and the strong-tunneling limit is always 
reached at sufficiently large energies. 



C. Interpretation of the quasiparticle statistics and 
calculation of the outgoing currents and shot noise 

In this subsection, we use the approach developed 
above to calculate the transport properties of the multi- 
point-contact junction in the limit of strong tunneling. 
Application of the voltage Vi to the Ith edge can be de- 
scribed simply as a shift of the incoming field (see, e.g., 
the Appendix) 



d in ) 



ki in ) 



nVit, 



(33) 



which can be understood as the evolution of the 
quantum-mechanical phase of the electron operator ipi 
(2) of the Ith edge. (In this work, we discuss only the dc 
transport properties of the junction biased by constant 
voltages Vi.) Numerical shift of the fields (pi does not 
affect their properties as quantum operators, and all the 
results for (pi discussed above remain valid in the presence 
of finite bias voltages. 

We consider first the situation with no backscatter- 
ing. Since the current ji in the Ith edge is related to 
(pi as ji(x,t) = —y f vidt(pi(x 1 t)/2ir, the outgoing currents 
are linear in the applied voltages and one can define the 
matrix of conductances through the current induced in 
the Ith edge by the voltage applied to the fcth edge: 
Gik = ji/Vk- The transformation of the fields (pi ac- 
cording to Eq. (32) means that the conductance matrix 
is: 



Gik — vo^JviVkS, 



Ik 



where o~q is the universal free-electron conductance e 2 /h 
equal to 1/2-7T in our units e = K = 1. Direct substitution 
of the matrices (31) into this expression gives 



Goo = v a Yl 



I/O + vi 



Cr 7 7 ^f 



J JJ ~ "3 



Gi 



V + Vj 



G tj = for j > i> 1, 



(34) 
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for 1 < j < % < n, 



2z/ 2A tt ^o - f; 



n 



^o + K? r- ^o + ^ 



where Gij denotes the conductance Gij in units of <jq. 
Charge conservation requires that ^ ; Gij = Vj, the con- 
dition that is indeed satisfied by Eqs. (34). 

Equations (34) describe a flow of charge along the 
edges in which the difference of currents in the edges 
and j is redistributed at the jth point contact ac- 
cording to its tunneling conductance <Jo(vj — Gjj) = 
2ctqVqVj j '(vq + Vj). Such a description of current flow 
in the system of several point contacts was first obtained 
in Ref. 10 under the assumption that quantum coherence 
of electron propagation along the common edge is sup- 
pressed, and electron distribution is equilibrated between 
the successive point contacts. The calculations presented 
above show that this description of the current flow re- 
mains valid in the quantum-coherent regime, when any 
external decoherence mechanisms are absent. 

We now turn to the regime of finite backscattering. As 
discussed above, each elementary process of backscatter- 
ing can be described as instanton tunneling generated 
by the operator ex exp{±i28j/\j} [see Eq. (26)]. If the 
tunneling occurs in the point contact j, the associated 
shift of the field <fij corresponds to the transfer of charge 
qj = 2/Aj = 2vqVj/{vq-\-Vj) (see, e.g., the Appendix). In 
the multi-point-contact junction, quasiparticle of charge 
qj produced in the Oth edge at the point Xj by instanton 
tunneling propagates along this edge and is split between 
edges i with i > j and the edge according to the conduc- 
tance matrix. Distribution of charges in backscattering 
processes that results from such a "current branching" 
at the point contacts enables us to interpret the quasi- 
particle Klein factors (27) of the dual model as exchange 
statistics of charges generated in the edge 0. Indeed, the 
absence of the energy gap at the edge of a FQHL makes 
it possible to produce arbitrary edge excitations. For an 
isolated edge 0, the exchange statistics of two excitations 
with charges p\ and p2 is set by its filling factor i/q to 
be "PxPiIvq. The two quasiparticles which tunnel in the 
nearest-neighbor contacts create charges qj and qj+i in 
the edge and the statistical angle for these quasiparti- 
cles, 2irqjqj + i/vQ is easily seen to indeed coincide with 
the statistical angle 2ttBjj + i determined by the Klein 
factors (27) through the statistical matrix By (22). If 
the quasiparticles tunnel at the point contacts that are 
not nearest neighbors, one needs to take into account 
that the charge created in the edge in the vicinity of 
the contact i by quasiparticle tunneling at the contact 
j < i is smaller than qj by the factor 



Yl [{vq - vi)/(y + n)] 



due to charge splitting at the intermediate point con- 
tacts. Equation (22) shows then that the matrix Bij can 
be written as Bij = rqiqj /vq, i.e., the statistics of quasi- 
particles can still be interpreted as the exchange statistics 
of charges in the common edge 0. 

This interpretation remains valid if some of the fill- 
ing factors v% are equal to each other. For instance, if 
the edges j and j + 1 have the same filling factor vq as 
the edge 0, the charge and statistics of the backscattered 
quasiparticles coincide with those of Laughlin's quasipar- 
ticles: qj = Bjj+i = vo, in agreement with previous 
results 3 . Our discussion here shows that in general the 
strong coupling between non-identical edges should pro- 
duce quasiparticles with both charge and statistics that 
are different from those of Laughlin's quasiparticles in 
the bulk. 

To find the backscattered currents quantitatively, we 
calculate first the charges injected into each edge by 
quasiparticle tunneling at the jth contact. Combining 
the two facts, that the fraction of the incoming current 
in the jth edge that goes out into the ith edge is equal to 
Gij/fj, and that the quasiparticle charge qj is created in 
the process of backscattering directly in the Oth edge, we 
see from Eqs. (34) that each backscattering event in the 
contact j creates the charge equal to Gij in the ith edge, 

i =^ j. This means that the backscattering current 7^ 
in the fth outgoing edge (I = ... n) is related to the rate 
of quasiparticle backscattering in the jth contact as: 
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To find the rates Jj , we need to calculate the dis- 
tribution of bias voltages for quasiparticles. In contrast 
to finding the bias (33) for electron tunneling, this task 
is not completely trivial, since in the strong tunneling 
limit, the bias voltage Vi applied to the Ith edge affects 
the quasiparticle bias Vj in all point contacts. The bias 
Vj can be obtained as the voltage-induced shift of the 
dual fields 6j in the arguments of exponents in the tun- 
neling Lagrangian (26). We find this shift from the linear 

relation between the incoming fields ^ and the fields 
9j (28). Extending the set of n fields 6j to the (n + 1)- 

component vector 8 = {<j>q ,81,82, ■■■, 8 n ) T we can write 
this linear relation similarly to Eq. (32): 

8 = S'<p {m) , S' = f { ^f {n -^ ■ ... -f {l \ (36) 



where the matrices T^> describe transformations of the 
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contacts. The part of T^ related to 

in P^ (30), while the part related to 6j follows directly 

from Eqs. (28) and (29): 
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(37) 



Since the voltages Vi shift the incoming fields accord- 
ing to Eq. (33), the transformation (36) shows that the 
voltage-induced shift of the quasiparticle fields 9j is 



6, 



Oj-^s'jiV^Vit. 



Obtaining the matrix elements S', t explicitly from 
Eqs. (37) and (31) we see that the voltage bias in the 
quasiparticle tunneling operator in Eq. (26) is: 
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(42) 



i.e., the charge is proportional to the fraction of the quasi- 
particle charge which reaches the Ith edge in the current 
branching process. The shot noise with the charge (42), 
generated at one point contact and fractionally split into 
another contact, is direct manifestation of coherent prop- 
agation of electrons between the point contacts. Equi- 
libration of electron distribution between the contacts 
would leave in each edge only the shot noise generated 
by the backscattering at the point contact formed by this 
edge. 



Inserting this result into the expression for the rate of 
quasiparticle tunneling from the Oth edge at the jth con- 
tact that follows from Eq. (26) we get: 
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(39) 



Finally, combining the strong-coupling branching of 
the currents according to the conductance matrix (34) 
with the backscattering contribution to the current, we 
see that the average outgoing current (/;) in the Ith edge 

is: 



(Ii, 



CO 



J2 G ki y i 



r (bsc)\ 



(40) 



In the limit of strong tunneling, the backscattering is 
weak and the average quasiparticle tunneling rate (/■ ) 
can be evaluated in the lowest order of perturbation the- 
ory in backscattering (26). For large quasiparticle bias 
voltages Vj one gets: 
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-sgnO^I^AD 
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This expression combined with Eqs. (35) and (40) de- 
scribes how the average outgoing current in the Zth edge 
approaches its linear large- voltage asymptotics. 

In this regime of weak backscattering, the quasiparti- 
cle tunneling produces regular shot noise, which at low 
temperatures is the only source of noise of the outgoing 
currents. 20 This mean that the correlators of the outgo- 
ing currents {!{) at zero frequency can be expressed as: 



(tt,4})U 



2]TAGyAG fc ,K/ 
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(41) 



where {...,...} denotes the anticommutator. In general, 
the current noise (41) is not proportional to the av- 
erage backscattered current, unless the contribution to 
backscattering of one point contact is dominant. If the 
contact m does dominate, the shot noise of the outgoing 
current in the Ith edge is characterized by the fractional 
charge equal to AG; m , 



V. TUNNELING BETWEEN MULTI-MODE 
FERMI-LIQUID RESERVOIR AND FQHL 

In this Section, we consider the situation with one spe- 
cific choice of the filling factors, Vj = 1, for j = 1 -=- n, 
in the general model of the multi-point-contact junction 
(Fig. 1). The multi-point-contact junction in this regime 
has been used before 9 ' 10 ' 12 to model an interface between 
the multi-mode Fermi-liquid reservoir and the FQHL. 
In the limit of strong tunneling, the reservoir-FQHL in- 
terface has all the features of a regular Ohmic contact: 
linear current-voltage characteristics and no shot noise, 
and establishes equilibrium between the reservoir and the 
edge of the FQHL when the number n of reservoir modes 
is large. Results of this section describe corrections to 
the Ohmic behavior caused by the weak quasiparticle 
backscattering which manifests itself through the shot 
noise and non- linear corrections to current. 

In the case of reservoir-FQHL interface, the v = 1 
edges of the general model of Fig. 1 represent one- 
dimensional free-electron scattering modes of the n-mode 
Fermi-liquid reservoir. When a voltage V is applied to 
the FQHL edge, the tunneling current It that flows be- 
tween the reservoir and the edge can be found as the 
sum of the outgoing currents [see Eqs. (35)and (40)] in 
the free-electron modes j — 1 -f- n: 



-aoAGooV - J2 G 03 lf(G j0 V) . 



(43) 
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The argument of the jth rate of backscattering I- here 
represents the quasiparticle bias voltage (38) in the jth 
contact. 

Perturbation theory in quasiparticle tunneling (26) of 
the dual model is justified when cither the temperature or 
all quasiparticle bias voltages are larger than the energy 
scale Tx — max J (T / F,/u) 1 ^ 1_2 / A ^ D of the crossover to 
strong tunneling in all point contacts. The condition on 
the bias voltages is always satisfied if the applied voltage 
V is sufficiently large: |V|[(1 - u )/{\ + u )] n > T X - 

Tunneling at large temperatures can be characterized 
by the linear conductance G of the reservoir-edge inter- 
face. We calculate this conductance by first obtaining 



the average rates (7- ) of backscattering in each contact 
in the first non-vanishing order of perturbation theory 
in tunneling (26) and then summing them according to 
Eq. (43). We find that with increasing temperature, con- 
ductance G approaches its saturation value G n , 

G n = a y [l-{q-l) n ], q =-^- = ^, (44) 

1 + f o -* 



reached in the absence of backscattering, as 
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Since q < 1, the first-order perturbation correction to 
G n (45) vanishes with increasing temperature as a power 
of temperature, G — G n ex 1/T 2 ^ q \ In the limit of 
large number of point contacts, n — > oo, this power-law 
correction is, however, suppressed by the same small fac- 
tor (1 — q) n which characterizes how G n approaches its 
equilibration value ctqi/ at n — > oo. One can show that 
correction to G n in the next non-vanishing order of the 
perturbation theory in backscattering is also suppressed 
by this factor. This suggests that G approaches the sat- 
uration value faster than any negative power of T. 

Following the same step as for the linear conductance 
(45) we obtain the average tunneling current in the large- 
voltage regime: 

(It) = G n V + { ~^t7 {V/D? {q - 1] (1 - <?) ( - 1)(2 *- 1) 
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For large n, the main contribution to the sum over j 
in this equation comes from the contacts at the end of 
the junction, where j ~ n. The overall magnitude of 
the average backscattered current (46) in the pertur- 
bative regime of weak backscattering is proportional to 
(1 — qf)( n_1 )( 29-1 ), and for uq > 1/3, does not vanish when 
n — > oo. If 1/q = 1/3, the exponent in this dependence is 
zero: 2q— 1 = 0, while for smaller v^ the exponent is neg- 
ative and makes the absolute value of the backscattered 
part of current (46) a growing function n. Equations (45) 
and(46) show also that the sign of the backscattering con- 
tribution to the tunnel current oscillates with the parity 
of the number n of point contacts, with backscattering 
increasing the tunnel current for even n. 

The shot noise of the tunneling current is given by 
the zcro-frcqucncy correlator ({It, It}} and similarly to 
Eq. (41) is generated by the backscattering part of the 
current (43): 



({7 T ,7 T })=^G 0i Go j ({4 0) ) /f}). 



(47) 



In the lowest non-vanishing order in the backscattering, 
the correlators on the right-hand-side of Eq. (47) are not 
equal to zero only for i = j. Non-zero diagonal corre- 
lators are given by the average backscattered currents, 
as in Eq. (41). In this case, summing the contributions 

2(7- ) to noise from each of the point contacts according 
to Eq. (47) we get 
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Similar to the average backscattered current, the main 
contribution to the shot noise comes from few contacts 
with j ~ n. If the tunneling amplitudes of these contacts 
are approximately equal so that we can neglect their vari- 



ations, W 2 ~ W 2 



the ratio between the noise and 



the average backscattered current is 



({It, It}) 
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We see that in a uniform junction with large number n 
of modes, the splitting of excitations at the point con- 
tacts complicates the relation between the backscattered 
charge and the shot noise. For instance, if v?, = 1/3, 
by a numerical coincidence Eq. (49) gives the noise-to- 
current ratio also equal to 1/3, while the backscattered 
charge is q = 1/2. For smaller vq, the ratio (49) does not 
correspond directly to cither q or vq. 



VI. CONCLUSION 

To summarize, this work provides a description of 
quasiparticle tunneling in the junction with geometry 
shown in Fig. 1, in which several point contacts are 
formed between single-mode FQHL edges with in gen- 
eral different filling factors Vj = l/odd and one common 
single-mode edge with filling factor i/q. This description 
extends to finite reflection in the contacts our previous 
theory 12 of strong electron tunneling in such a junction. 
New quasiparticle physics emerging form the multi-point- 
contact geometry is coherent splitting of quasiparticles 
at all point contacts downstream the common edge v$ 
from the point contact where the quasiparticles have been 
created. Quasiparticles are split at each contact in pre- 
cisely the same way as the incoming edge currents in the 
strong tunneling limit. This means that the branching 
ratios of the edge currents determine the quasiparticle 
charge and statistics. The statistics is given by the per- 
mutation relations between the backscattering operators 
producing quasiparticles and makes the whole process 
of coherent splitting at successive point contacts causal. 
The quasiparticle charges determine the magnitude of the 
shot noise of the outgoing current in each edge. In the 



non-uniform junctions, in which one point contact domi- 
nates the backscattering, the shot noise in any given edge 
is directly related to the charge split in this edge in the 
branching process - see Eq. (42). In uniform junctions, 
the charge-noise relation is less direct [see Eq. (49) and 
the example discussed next to it] because of averaging 
over different point contacts. 

In the case of a junction between the Fermi-liquid 
reservoir and the FQHL edge, Vj = 1 for j = 1 -=- n, 
the quasiparticle backscattering gives corrections to the 
reservoir-edge equilibration reached with increasing num- 
ber of reservoir modes, n — > oo, in the absence of 
backscattering. For temperatures T much larger than 
the bias voltage V, these corrections are uniformly sup- 
pressed by the factor \(v — l)/{v + 1)]", exponential in 
the number n of modes. In the opposite regime, T«y, 
the corrections, i.e. the non-linear junction conductance 
and current shot noise, do not vanish, however, even for 
n — ► oo. 

Description of the charge and statistics of the edge ex- 
citations in the multi-point-contact junction, which is de- 
veloped in this work for simple single-mode edges in the 
junction geometry with one common edge for all point 
contacts, should also be important for more complex 
multi-mode edges and more general junction geometries. 
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APPENDIX A: ONE POINT CONTACT 

In this Appendix, we review the known results 8 ' 18 for 
strong tunneling in one point contact between two edges 
with different filling factors using the approach similar to 
that employed in the main text. For consistency of nota- 
tions, we number the two edges and j, with the filling 
factors vq and Uj . The Hamiltonian of one point contact 
can be written as in the model described in Sees. II and 
III 



n = 



dx (£-[(d^ ) 2 + {dxfaf} + [Vo Po + V jPi 



(Al) 



The two differences with Sees. II and III are that now we 
explicitly include the bias voltages V~o and Vj applied to 
the edges, and omit the Klein factors in the tunnel part 
of the Hamiltonian, since they are irrelevant in the case 
of one point contact. (The Klein factors £o£j and £i£o in 
the forward and backward tunneling terms commute with 



each other and are constants of motion. This means that 
they can at most produce an irrelevant constant shift of 
the fields 4>.) 

The tunneling field ipj is introduced according to 
Eqs. (29) which can be viewed as the rotation in the 
space of the fields <\>: 
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cos dj , — sin dj 
sin {>j , cos -dj 



by the angle #j defined through the relations: 
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Performing this rotation in the Hamiltonian (Al), 
in particular rotating the "vector" of bias voltages 
^Jvq Vo , J~v]Vj , we see that the Hamiltonian Tix of the 
tunneling mode Lpj can be separated from other terms: 

-^\^{dxVj? + —dx^j] - Uj cos[Xjipj(x = 0)] , 



V=V - V, 



(A2) 



From now on we omit the index j of all variables. Hciscn- 
berg equation of motion for if obtained from the Hamil- 
tonian (A2) and the commutation relations (5) is: 

V(t) 

(d t + vd x )f(x,t) = h u7rAsin[A<p(0,t)]sgn(a;) . 



Separating the bias voltage: 

1 /■* 
(fi(x,t) = tpo(x,t) - - / dt'V(t') 

we can write the solution of Eq. (A3) as 
2irm 



(A3) 
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6(t-x/v) + n\Q(t), x <0, 
X(t - x/v) - n\Q(t), x > , 



(A5) 



where an arbitrary constant 2irm/\ is written in this 
form for later convenience, and Q(t) is the operator of 
charge transferred through the point contact: 



Q(t) = -usm[\ip(0,t)]. 



(A6) 



The function 9(t) in Eq. (A5) is an arbitrary operator 
which has the meaning of fluctuations of the field (f(x, t) 
incident on the point contact, and x will be determined 
from the final solution of Eq. (A3). 

Equation (A3) implies that (p(x, t) is continuous as a 
function of x at x — 0, i.e., 



<A)(0,i)- 
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= 6{t) + nXQ(t) = X (t) - n\Q{t) . (A7) 
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In the strong tunneling limit u — > oo, we assume that the 
value of <p(0, t) is very close to one of the minima 2irm/\, 
m = 0,±1, ... of the tunneling term in the Hamiltonian 
(A2) and we can linearize the sine in Eq. (A6). As we 
will sec in the end of the calculation, this assumption is 
indeed correct. Applying it to the first of the Eqs. (A7), 
we get the following simple equation for ipo(Q,t): 



that the large-scale distribution of the field <f(x, t) cre- 
ated by such a process of instanton tunneling occurring 
at time ti coincides essentially, as it should, with the re- 
tarded Green's function discussed in Sec. II: 



, n 2"" / 

y{x,t) = -— 0(i 



ii)sgn[x - v(t - ti)] 



(A13) 



<po{0, t) = -uTT\[\tp (0, t)- I dt'V(t')\ + 9 . (A8) 



Solution of this equation under the condition 9, V <C u 
appropriate for the strong-tunneling limit gives 



<p(o,t) 
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We see that up to small terms of the order of 1/w, 
</?(0, t) is indeed pinned to one of the minima 27rm/A of 
the tunneling energy. Inserting Eq. (A9) into Eq. (A6) 
we see that the tunnel current in the point contact is: 



OM.^vw-M 



(A10) 



The first term on the right-hand-side of this equation 
gives the strong-tunneling conductance of the point con- 
tact, while the second term represents the tunnel cur- 
rent induced by the incident field fluctuations. Finally, 
combining Eq. (A9) and (A10) with the continuity con- 
dition (A7) we can determine the transmitted field %(£) 
in Eq. (A5): 



x(t) 



A 



dt'V(t') - 9(t) . 



(All) 



(Equations (A10) and (All) include only the terms that 
do not vanish in the limit u — > oo.) Comparison of 
Eq. (All) with Eqs. (A4) and (A5) shows explicitly that 
both the voltage contribution to <p(x, t) and incident fluc- 
tuations 9{t) change sign at the point contact in the 
strong-tunneling limit. This means that the field 6j(x) 
defined by Eq. (28) of the main text indeed propagates 
freely in the strong tunneling limit. 

The last remark is that Eq. (A5) enables us to ob- 
tain directly the charge of backscattered quasiparticles. 
Indeed, instanton tunneling that describes the backscat- 
tcring changes the number m of the minimum the field 
f(x, t) is pinned to, affecting in the process the charge 
transfer through the point contact. Equation (A5) shows 
that the transition from one tunneling minimum to an- 
other changing m by 1 happens if the transferred charge 
changes by the amount 



SQ = - 2 



(A12) 



which coincides with the junction tunnel conductance in 
units of <jq. One can also see from Eqs. (A5) and (A7) 
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